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DYNAMICS  OF  LIQUID  PROPELLANT  GUN  LOADING  SYSTEMS 
Annual  Progress  Report,  FY75 


I.   INTRODUCTION  AND  BACKGROUND 

A.  LIQUID  PROPELLANTS  FOR  ADVANCED  NAVY  GUN  SYSTEMS 

With  continued  trends  toward  more  sophisticated  weapons  platforms 
to  meet  the  projected  threats  of  a  hostile  environment,  improved  gun 
systems  will  be  required.  These  systems  of  the  future  will  have  to 
have  a  total  weight  less  than  existing  gun  systems  of  the  same  cali- 
ber, a  less  complicated  logistics  train,  improved  capabilities,  and 
they  will  -have  to  be  cost  effective.  Various  analytical  and  experi- 
mental investigations  conducted  over  a  period  of  approximately  thirty 
years  have  demonstrated  that  the  potential  advantages  of  liquid  propel 
lant  guns  (LPG's)  are  sufficient  to  warrant  their  serious  considera- 
tion as  candidate  elements  in  advanced  naval  gun  systems. 

The  following  remarks  are  intended  to  illustrate  some  of  these 
potential  advantages.  By  eliminating  the  cartridge  case,  both  air- 
borne and  surface  platforms  will  be  able  to  carry  more  shots  per 
pound  and/or  a  more  varied  assortment  of  weapons.  At  the  present 
time,  shipboard  guns  using  semi -fixed  ammunition  require  a  separate 
magazine  for  powder  in  relatively  close  proximity  to  the  gun.  By 
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using  a  liquid  propel lant,  the  propellant  tanks  could  be  located 
anywhere  in  the  ship  without  great  difficulty,  thus  allowing  the 
ship  designer  more  flexibility.  With  a  standardized  propellant, 
different  caliber  guns,  which  now  use  both  semi-fixed  and  fixed 
ammunition,  could  be  supplied  from  the  same  tank,  eliminating  the 
separate  spaces  needed  to  conform  with  ordnance  storage  compatibility 
requirements.  Redistribution  of  propellants  onboard  could  be  accom- 
plished by  means  of  pumping  rather  than  manual  labor  or  an  elaborate 
mechanical  transfer  system.  Propellant  retrograde  materials,  e.g., 
pallets,  powder  cases,  and  powder  cans,  would  be  eliminated.  Under- 
way replenishment  and  rearming  evolutions  would  consist  of  taking  on 
projectiles  and  liquid,  thus  significantly  reducing  time,  cost,  and 
effort.  Ship  stability  might  be  improved  by  the  elimination  of  propel 
lant  weight  in  ready-service  spaces.  A  ship  sustaining  battle  damage 
could  easily  improve  her  stability  by  redistributing  propellant  or 
pumping  it  over  the  side.  Increased  firing  rates  may  be  attainable. 
It  is  because  of  this  prospectus  and  improved  technology  in  the  field 
of  monopropellants  and  bipropellants  for  gun  applications  that  a  re- 
search and  development  effort  has  again  been  initiated  in  the  liquid 
propellant  gun  (LPG)  area  of  ordnance  development. 

B.  REVIEW  OF  PREVIOUS  DEVELOPMENTS 

In  order  to  establish  a  point  of  departure  for  current  develop- 
mental efforts  in  LPG  technology,  it  is  useful  to  briefly  review  the 
experiences  documented  by  previous  investigators.  A  review  of  LPG 
history  was  conducted,  but  the  following  remarks  refer  only  to  un- 
classified literature. 


Interest  in  LPG  systems  in  the  United  States  was  evident  as  early 
as  1945  [1]  when  the  Army  Office  of  the  Chief  of  Ordnance  conducted 
some  preliminary  calculations  of  the  performance  expected  from  liquid 
propel lants.  It  was  concluded  that  artillery  performance  could  be 
greatly  improved  by  using  liquid  propellants  instead  of  solid  pro- 
pellants.  Initial  experimental  investigation  was  in  the  area  of  bi- 
propellants  and  feasibility  demonstrations  were  conducted  with  a 
smooth  bore,  caliber  .50  gun.  These  early  tests  resulted  in  surpris- 
ingly high  muzzle  velocities  which  seemed  to  spur  increasing  interest. 
It  was  discovered  that  to  inject  the  propel lant  properly,  an  injection 
pressure  in  excess  of  chamber  pressure  was  required  [1J.  To  circum- 
vent this  problem,  pre-loaded  cased  rounds  were  used.  Large  variations 
in  the  ballistic  behavior  without  repeatability  were  noted  [2]. 

In  1948,  investigation  into  the  feasibility  of  regenerative  injec- 
tion of  bipropel lants  was  begun  [3].  The  regenerative  injection  tech- 
nique utilizes  a  differential  area  piston  to  extract  energy  from  the 
burning  propellant  to  power  the  injection  process.  It  was  determined 
that  regenerative  injection  was  feasible.  However,  continued  study 
through  1950  with  caliber  .50  guns  emphasized  the  need  for  high  pres- 
sure static  and  dynamic  seal  improvement.  It  was  also  noted  that 
lower  mixing  intensities  gave  lower  muzzle  velocities  and  excessive 
round-to-round  variations  were  encountered  with  dilute  bipropel lants 
[4].  On  the  positive  side,  liquid  propellants  produced  less  bore  wear 
than  the  existing  solid  propellants. 

Concurrently,  work  was  started  on  a  37-mm  gun,  based  on  the  caliber 
.50  work  [5,  6].  The  37-mm  gun  had  regenerative  injection  with 


nitrogen  pressure  to  start  the  process.  Feasibility  of  using  the 
regenerative  process  in  37mm  guns  was  proven;  however,  further  dev- 
elopment in  the  area  of  seals  was  again  indicated. 

The  Navy  became  interested  and  began  investigation  of  LPG  systems 
in  1949  with  an  emphasis  on  improved  performance  for  aircraft  guns 
[7,  8].  Larger  guns  with  higher  velocities  at  the  same  size  and 
weight  were  also  desired.  Naval  Ammunition  Depot,  Seal  Beach,  Calif- 
ornia, started  research  and  development  of  a  40mm  monopropellant  gun. 
The  monopropellant  developed  seemed  to  satisfy  all  of  the  desired 
propellant  criteria,  however,  for  their  initial  gun  tests,  a  bipropel- 
lant  was  used  in  a  pre-loaded  casing.  At  the  same  time,  separate 
research  was  being  done  in  the  area  of  propellants  themselves  [9]. 

It  should  be  noted  that,  as  of  1950,  the  majority  of  the  inter- 
est in  liquid  propellant  guns  was  in  the  Army.  They  were  mainly  in- 
terested in  LPG  technology  vice  development,  with  a  major  interest  in 
the  end  application  to  antiaircraft  weapons.  The  high  muzzle  veloci- 
ties and  high  rates  of  fire  anticipated  could  readily  satisfy  the 
requirements  for  this  role.  Another  use  contemplated  for  LPG's  was 
aircraft  machine  guns,  because  of  the  high  rates  of  fire  and  the 
elimination  of  cartridge  cases. 

With  the  commencement  of  hostilities  in  Korea,  the  Army's  LPG 
interest  shifted  to  the  area  of  infantry  and  tank  weapons.  The  hyper- 
velocities  and  greatly  increased  number  of  rounds  per  reload  were 
extremely  desirable  characteristics  for  infantry  weapons.  A  program 
specifically  designed  to  employ  liquid  propellants  in  an  existing 
90mm  gun  as  a  tank  weapon  was  started  in  1950  [10].  This  retrofit 


program  was  aimed  at  improving  the  armor  penetration  of  Army  guns  in 
Korea.  The  concept  was  to  fit  into  a  standard  cartridge  case  the 
components  of  a  hypergolic  bipropellant  system  necessary  to  achieve 
muzzle  velocities  in  excess  of  a  standard  round.  Hypergolic  bipropel- 
lants  are  designed  to  ignite  spontaneously  upon  mixing.  For  a  liquid 
propel! ant  gun,  this  implies  that  no  initiator  is  needed  for  ignition. 
These  hypergolic  bipropellants  require  an  injection  or  mixing  process. 
In  the  event  of  accidental  contact  of  fuel  and  oxidizer,  a  fire  would 
inevitably  occur.  Separate  fuel  and  oxidizer  containers  were  placed 
in  the  cartridge  case  with  a  mechanism  to  effect  mixing.  After  two 
years  of  work,  the  design  objectives  had  not  been  reached;  only  a  half 
charge  of  propel lant  could  be  used  without  over-pressurizing  the  gun. 

Actual  small -caliber  design  studies  of  caseless  guns  were  started 
in  the  early  1950's.  The  major  emphasis  in  this  field  was  on  aircraft 
machine  guns.  A  caliber  .60  gun  was  designed  in  1953  incorporating  a 
regenerative  injector  monopiston  to  introduce  both  fuel  and  oxidizer 
in  equal  amounts  into  the  chamber  [11].  The  design  rate  of  fire  was 
2,000  rounds  per  minutes.  However,  the  projectile  and  propellant  feed 
systems  interacted  in  such  a  way  that  projectile  feed  could  inhibit 
propellant  loading.   It  was  recommended  that  a  revolver  type  feed 
system  be  investigated.  The  reproducibility  of  ballistic  character- 
istics, misfire  procedures,  overpressurization  of  the  gun,  propellant 
ignition,  and  ullage  problems  continually  plagued  these  early  efforts 
in  liquid  propellant  gun  research. 

Until  1955,  LPG  systems  consisted  primarily  of  two  general  classes, 
preloaded  and  injection  guns.  Bipropellants  and  monopropellants  have 


been  employed  in  each  system  with  varying  degrees  of  success.  Pre- 
loaded guns,  though  simple  in  design,  were  extremely  sensitive  to 
loading  density  and  ignition  variations.  Injection  guns,  though 
relatively  complex,  offered  a  positive  means  of  combustion  control. 
In  an  attempt  to  utilize  the  advantages  of  the  various  systems  pre- 
viously investigated,  a  project  was  started  in  1955  to  combine  the 
two  systems  into  a  single  monopropellant  system  approaching  the  sim- 
plicity of  a  preloaded  gun,  yet  having  the  advantage  of  combustion 
control  through  injection  [12].  A  40mm  model  was  designed,  featuring 
regenerative  injection.  A  low-pressure  static,  high-pressure  dynamic 
seal  configuration  was  adopted  to  circumvent  earlier  problems  encountered. 
Marginal  performance  was  obtained  with  the  system. 

At  the  same  time,  the  Navy  funded  research  in  the  area  of  monopropel- 
lant and  bipropellant  injector  systems  for  a  30mm  gun  for  aircraft  appli- 
cations. This  design  incorporated  a  breech  bleed  feature  to  reduce  ullage, 

A  40mm  development  program  was  started  as  a  prototype  for  a  3  to  5 
inch  rapid  fire  gun  [13].  The  40mm  was  selected  for  testing  purposes  as 
a  caseless  gun.  This  effort  was  apparently  terminated  after  the  gun  was 
damaged. 

A  number  of  projects  were  in  progress  in  the  mid-50' s  which  were 
aimed  at  producing  a  90mm  tank  gun  [14,  15,  16,  17].  These  programs 
continued  to  strive  for  hypervelocities  by  employing  regeneratively 
injected  bi propel! ants.  An  impressive  increase  in  gun-tube  life  was 
anticipated;  however,  the  injection  system  was  somewhat  large  for  tank 
application.  At  this  point  (1957),  development  aimed  at  major-caliber 
end  items  for  Army  use  appears  to  have  stopped,  although  some  small -arms 


studies  continued.  Among  the  reasons  for  this  curtailment  of  research 
would  appear  to  be  the  movement  to  a  peace-time  economy  with  the  cessa- 
tion of  hostilities  in  Korea  and  the  growth  of  missile  technology. 

In  1968,  programs  were  started  by  the  Navy  which  addressed  the 
development  of  a  safer  aircraft  machine  gun  with  improved  performance 
[18].  The  use  of  non-hypergol ic ,  preloaded,  caseless  liquid  bipropel- 
lants  was  the  course  of  action  selected.  After  encouraging  initial 
results,  the  program  was  continued  with  emphasis  placed  on  developing 
a  20mm  burst-fire  gun.  These  programs  utilized  red  fuming  nitric  acid 
(RFNA)  as  a  propel lant  ingredient.  One  of  these  programs  investigated 
a  rotating  four-barreled,  cam  operated  system  designed  for  a  rate  of 
fire  of  1000  rounds  per  minute.   It  employed  a  pressurized  injector 
feed  system  having  an  injection  period  of  11.3  milliseconds  [19].  A 
computer  simulation  indicated  the  system  would  perform  satisfactorily; 
however,  a  working  model  was  not  constructed.  A  similar  program  con- 
ducted simultaneously  also  employed  a  multi-barreled  design  [20].  A 
single  barrel  test  gun  was  manufactured  and  fired  in  single  shot  tests, 
In  the  multi-barreled  design,  the  projectile  would  be  seated  by  propel - 
lant  pressure.  A  nitrogen  pressurized  accumulator  system  for  the  prop- 
ellant  was  employed  to  maintain  propellant  pressure.  For  the  inert 
tests,  water  was  used  as  the  propellant  simulator. 

C.   SCOPE  AND  OBJECTIVES  OF  THE  PRESENT  STUDY 

The  overall  goal  of  this  investigation  is  to  both  analytically  and 
experimentally  simulate  and  study  the  behavior  of  a  liquid  under  con- 
ditions analogous  to  those  existing  in  a  rapid-fire  LPG  feed  system. 
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The  results  of  the  study  are  intended  to  be  useful  in  the  identifica- 
tion of  upper  and  lower  bounds  for  such  LPG  design  and  performance 
parameters  as  time-to-load,  injection  supply  pressure,  injection  sys- 
tem configuration  (line  lengths,  diameters,  accelerating  masses,  and 
viscous  and  minor  losses),  ullage,  charge-to-mass  ratio,  caliber  size, 
and  projectile  mass.   In  addition,  the  investigation  has  been  designed 
to  detect  and  characterize  "anomalous"  behaviors  such  as  chamber  pres- 
sure oscillations  during  loading,  existence  and  consequences  of  cavita- 
tion in  the  liquid  slug,  and  various  sources  of  delay  in  the  loading 
system.  It  has  been  felt  that  laboratory  observation  of  effects  such 
as  these  will  prove  to  be  essential  in  understanding  the  behavior  of 
LPG  prototype  systems  and  future  operational  devices. 

D.  APPROACH  AND  PROGRESS 

The  investigation  has  taken  the  classical  approach  of  analysis  com- 
plemented by  experimental  tests.  The  analytical  model  has  been  con- 
structed with  a  view  towards  simplicity  while  maintaining  flexibility 
to  refine  the  analysis  as  required  to  meet  design  requirements  and  to 
incorporate  new  information  as  LPG  technology  advances.  The  analytical 
model  in  its  current  form  assumes  one-dimensional  unsteady  flow  of  an 
incompressible  propel lant.  This  latter  assumption  is  viewed  as  the 
most  critical  deficiency  in  the  analysis  and  is  a  recommended  area 
for  follow-on  efforts.   It  will  be  seen,  however,  that  the  analytical 
model  in  its  present  form  is  adequate  for  use  in  obtaining  prelimin- 
ary design  guidance. 

The  experimental  portion  of  the  study  has  been  primarily  designed 
to  support  and  validate  (hopefully)  the  analytical  model.  Accordingly, 
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no  special  effort  has  been  made  to  emulate  existing  or  envisioned  LPG 
designs.  Off-the-shelf  hardware  has  been  used  in  the  construction  of 
the  laboratory  LPG  simulator  so  that  the  performance  of  the  device  is 
somewhat  less  than  might  be  obtained  from  an  optimally  designed  system. 
The  analytical  model  accurately  predicts  this  "off-design"  performance 
and  points  the  way  toward  the  areas  of  redesign  that  are  likely  to  be 
most  beneficial.  This  redesign  is  another  area  recommended  for  the 
expenditure  of  future  efforts. 

Progress  to  date  may  be  summarized  as  follows: 

1.  An  analytical  model  of  LPG  loading  systems  has  been  constructed 
and  exercised.  The  model  has  been  shown  to  adequately  predict  the 
overall  performance  of  LPG  loading  systems  (time-to-load)  and  the 
dominant  factors  to  be  considered  in  the  design  of  LPG  systems. 

2.  An  LPG  laboratory  simulator  has  been  designed  and  built  and 
has  been  shown  to  be  a  useful  tool  for  the  evaluation  of  analytical 
predictions  and  for  the  identification  of  unpredicted  loading  system 
behavior. 

3.  The  analytical  model,  verified  by  the  experiments,  has  been 
used  to  provide  design  guidelines  for  the  estimation  of  LPG  loading 
system  performance  as  a  function  of  various  design  parameters. 

These  developments  are  discussed  in  detail  in  subsequent  sections 
of  this  report.  The  report  is  concluded  with  a  summary  of  results  and 
recommendations  for  future  work  in  the  subject  area  of  investigation. 
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II 


ANALYTICAL  MODEL 


Analytical  efforts  have  been  aimed  at  developing  a  model  for  the 
behavior  of  a  liquid  as  it  is  pumped  from  an  injector,  through  a  system 
of  connecting  lines  and  valves,  into  the  breech.  Although  continuing 
analytical  work  will  endeavor  to  treat  this  behavior  in  some  detail  (in- 
cluding, for  example,  propagation  of  pressure  waves  and  local  occurrences 
of  cavitation),  the  complexity  of  the  process  has  dictated  a  quasi-lumped 
parameter  approach  for  initial  studies.  As  shown  in  the  sketch  below, 
the  fluid  is  considered  to  be  made  up  of  three  slugs  occupying  the  in- 
jector, connecting  lines,  and  breech.  The  present  dynamic  analysis  takes 
into  account  the  transfer  of  fluid  from  one  region  to  another  in  response 
to  a  step  input  in  pressure  at  the  injector  head  end. 
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In  order  to  conform  with  current  design  concepts,  including  the  design 
of  the  laboratory  simulator,  the  system  is  modeled  in  the  "projectile- 
ram  mode"  in  which  the  breech  is  the  same  diameter  (nominally)  as  the 
bore.  The  model  may  be  easily  extended  to  the  "flow-thru"  mode  of 
chamber  filling  by  (in  the  first  approximation)  simply  neglecting  the 
mass  of  the  slug. 

A.  CONSTITUENT  RELATIONSHIPS 

1.  Injector  inertia:  The  loss  in  pressure  across  the  moving 
injector  mass  (see  sketch)  may  be  obtained  from  the  force  balance: 


PA  -  P.A.  =  m  x.+  k-x.  +  k  x.  +  F 
r  r         l  i    pi   f  p  i    sp  l    p 


r 

D     P  -i 
r     r  _^ 


n 


~P.   D. 
l    i 


^__t 


(1) 


2.  Losses  in  the  liquid:   In  this  region  there  must  be  a 

pressure  drop  to  overcome  both  the  inertia  and  drag  of  the  liquid, 

Thus,  in  the  injector 

U,  -  x,)  x/ 

(2) 


(L.  -  x.)  x.2 

APi  =  ^Li  -xi)*i  +  pfi~n 


D 


in  the  connecting  line 

L  x  2 

A?t   =  PLA  +  Pf£  ^  -\ 

and  in  the  chamber 


APC  =  pXsXs   +  pf  Ji  !| 

c 


(3) 


14 


In  addition,  minor  losses  can  be  taken  into  account  by  writing 

•  2 

where  k. 


mi    •  2 
2 
For  example,  for  the  drop  due  to  the  sudden  contraction  from  the 

injector  to  the  connecting  line,  we  have  (i  =  1) 

AP     =  Ifi..   2  =  3^  2   (XA 
Ahl       2  2  xl         2  2  xs      I  •    / 

vxs' 

andkml  =  f(^)     * 

Likewise,  for  orifices  occurring  in  the  connecting  line,  we  would 
obtain  (i  =  2),  for  each  such  occasion, 

2  n         /  x2   •    v2 


2      irX"        2  xl         \  CX     /l  •    /    2  xs 


do/  \dc/\x 


andk^  =  fe)2(r)2- 

s 

Invoking  the  incompressible  assumption  the  velocities  can  all 

be  written  in  terms  of  x  .  Thus,  since  the  volume  flow  from  the 

injector  must  equal  the  inflow  to  the  chamber, 

A.  x.  =  nAJL  =  A  x  (5) 

i  i     I  I        c  s  K 

whence                    A.,  x\   =  nA£x£  =  Acx's  (6) 

and                          A.x.   =  A  x  (7) 

lies  K    ' 
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Combining  Eqs.  (2)  through  (7),  the  total  pressure  drop  across  the 
liquid  may  be  written, 


(8) 


px. 


f. 
i 

D. 
l 


t 


Pressure  drop  from  reservoir  to  r.lug: 
By  simple  manipulation,  we  have 


P   -  P'=  (PA   -P. A.)  Tr-        +  (P.  -  ?')   +  P   (I-tt1) 

r    s   K   r  r         i  l '  A.     v  i    s y    r  \    A.  / 

so  that  after  applying  Eqs.  (5)  through  (7)  to  Eq.  (1),  combining 
with  Eq.  (8),  and  collecting  terms,  the  total  drop  from  reservoir 
to  slug  is 


p -p: 


-  =  1    -  Ps  =  k1   +  k2  x   (H]   +  H^)  +  k3  x^   (H3  +  H4x) 


+  k4  x  +  k5x  +  k6 


p: 


where  Ps  =  ~  ,  x  =  -^  ,  and 
r  sm 


k1  =  i 


ft) 


m     x 
k     _  _i sm 

K2  "  P     A 
r     c 

t     m     x     2 
v     =  1      s  Asm 
K3  "   2     PAD 
r  c  c 


(t) 


=  JlL 


PrA 


=  _§JL 


x         /  D    \ 
sm    [    _c  ) 

c        V  Di  / 


P     A 
r     c 


kr  = 


P  A 
r  c 


(fe) 


(9) 
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c_ 

m 


px   A 

(the  charge- to-mass  ratio) 


m 


D 


4 


ms     \Di/        ms 


1  - 


D. 
l 


L. 

l 


sm 


m, 


Ti        x         1  D.  /  2       x        I 

sm     \    i  /         n  sm   \ 


£ 


D 

+  —  5X- 

x       vmi 
sm    i 


■V',(!f 


4.  Drop  across  the  slug : 

Again,  writing  a  force  balance, 

P'A     =  m     x     +  k.    x     +k     x     +F     +  P.    A 
sc         ss         fss         sss         s         be 


or     P„  =  k0  x  +  k-,  x  +  k0  x  +  kn 
s         2  7  8  9 


where 


(10) 


k,    x  k    x  F  p. 

i,  ts     sm         ,      _     ss  sm       v  s       ,  __t^ 

K7         P^  A„       '   K8  '     P^  kr   '   N9       PA         P^   * 
re  re  r  c         r 

Equations  (9)  and  (10)  are  the  governing  differential  equations 

for  the  slug  position  and  slug-face  pressure  as  a  function  of  time. 

In  subsequent  sections,  we  have  developed  solutions  for  these 

equations  under  the  initial  conditions 

x(0)  =  x(0)  =  0 
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B.   SOLUTIONS  TO  THE  GOVERNING  EQUATIONS 

Equations  (9)  and  (10)  have  been  combined  and  programmed  for 
solution  on  the  analog  computer.  This  mode  of  solution  has  been 
chosen  because  of  the  ease  of  parameter  variation  made  possible  by 
the  analog  method.  The  analog  computer  solution  is  outlined  in 
Section  IV  of  this  report.  In  this  subsection,  several  closed  form 
solutions  are  presented  as  they  are  especially  useful  for  the  spec- 
ial cases  treated. 

1.  Solution  of  the  semi-linear  case.   Combining  Eqs.  (9)  and 
(10),  we  obtain,  upon  elimination  of  P, 
yO+H-j+^x)  +  k3x2(H3+H4x)  +(k4+ky)x  +  (k5+k8)x  =  1  -  (k-, +k6+kg)   (11) 

This  expression  is  non-linear  due  to  the  fluid  friction  term 
(preceded  by  k3)  and  the  terms  accounting  for  the  shift  of  fluid 
from  the  injector  to  the  chamber  (terms  preceded  by  H2  and  FL).  An 
important  simplification  can  be  made  if  the  system  design  parameters 
are  such  that 

H4        H2 

and  if  kr  +  k~  <<  1 . 

The  first  set  of  inequalities  requires  that  the  difference  between 
the  velocity  and  acceleration  of  the  fluid  in  the  chamber  and  of  that 
in  the  injector  is  small  or  that  piston  mass  and  connecting  line  ef- 
fects dominate  the  system.  The  second  inequality  neglects  spring 
effects  in  the  piston  and  slug  behavior.  Both  of  these  conditions 
represent  real  (but  not  necessarily  desirable)  situations  and,  in 
fact,  are  close  approximations  to  the  conditions  of  the  experiments 
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described  in  Section  III  of  this  report.  With  these  simplifications 
we  have, 

k2(l  +  H^x  +  k3  H3  x2  +  (k4+k7)x  =  1  -  (k^kg+kg) 

x(0)  =  x(0)  =  0 
This  is  a  Riccati  differential  equation  with  the  solution 


(13) 


x  =  A  In   cosh  B(t+to)  -  Ct 
cosh  Bt 


x  =  AB  tanh  B(t+t  )-  C 
o 

x  =  AB2  sech2B(t+t  ) 

o 

k2(1+Hl}  /T~  k4+k7 

WHere  A  =  ~W~  '   B  =     2k2(l+Hl)    '   C  =  2k^3 

q  =  (k4+k7)2  +  4  [l  -  (k^kg+kg)]  (k3H3) 
1  -1    (k4+k7) 

Vftanh      "7T 

The  variation  of  pressure  P  with  time  may  be  obtained  by  substituting 

Eq.  (13)  into  Eq.  (9)  or  Eq.  (10),  but  this  rather  laborious  procedure 

has  been  omitted  in  deference  to  the  slavishness  of  the  computer. 

A  most  useful  application  of  Eq.  (13)  is  to  use  it  to  estimate 

the  time-to-load,  t*.  When  t  =  tr,  x  =  x       or  x  =  1 ,  so  that 

f  f   s   sm 

cosh  B(tf  +  t  ) 

1  =  A  In     t-dt ~  -  Ct. 

cosh  Bt         f 
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Under  the  further  approximation  (also  applicable)  that  piston  and 
slug  damping  forces  are  negligible,  we  have  C+o  and  t->o,  giving 


where 


B1  = 


t,  =   g-,  cosh"  exp  -n- 


|[1  -  (^  +  k6  +  kg)]  (k3H3)  } 


(14) 


V2 


2k2(l  +  H^ 


Each  of  the  terms  in  this  expression  contain  design  parameters.  In 
particular,  it  can  be  immediately  seen  that 


tf  a  B, 


'2  m  D 
s  c 

P  A 
,     r  c 


1/2 


1   +  H 


r 


{(thu 


fV 


|V2 


(15) 


where  H-,  and  H3  are  defined  with  Eq.  (9)  and  P.  is  an  effective 
back  pressure  due  to  static  drag  and  back  pressure  given  by 

F    F 

P  =  — 2-  +  —  +  P 
Pd   A,   Ac   Pb 

The  value  of  B'tf  given  by  Eq.  (14)  is  plotted  in  Fig.  1.  These 
expressions  will  receive  additional  attention  in  the  discussion  of 
the  experiments. 
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Btf 


0 


0 

• 

2 

Figure    1. 


tf 


|_2k7TT 


+   Hx) 


o 

1.5 

1    +         H2 

1    +    H- 

Figure    2. 
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2.  Other  closed  form  solutions:  Other  less  realistic  simplifica- 
tions can  be  made  to  Eq.  (11)  and  these  are  listed  here  for  completeness 
Further  details  are  given  in  Appendix  B. 

a.  Inertia  effects  only:  Here, 


where  a  =  1 


k2  x  (1  +  H-j  +  H2x)  =  a 


2 

<ki  +  ke +  k9»  -  (bt)    -  rv 


The  solution,  developed  in  Appendix  B,  is 

H, 
t  =    f  1   + 


1/2 


(]  +T^TX) 


2k2   (1   +  H^/l   +  H^ 


D(u) 


where  D(u)   is  Dawson's  integral   [21] 


D(u)  =  e 


jr. 


e      dt 


and,  u  =     [lr\  (1   +  ]   +  ^    x)] 


1/2 


Setting  x  =  1  we  have,  for  the  ram  time 


2k2(l  +  H^ 


1/2 


1  +  H 


=  1  + 


1/2 


1  +  H1 /\  H2 


D(u)   (16) 


This  result  is  plotted  in  Fig.  2  for  representative  values  of 

Ho 
1  + 


1  +  H]  * 
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Note  that  H  is  a  measure  of  the  nonlinearity  contributed  by  the 
difference  in  acceleration  of  the  liquid  between  injector  and  piston 
(Hp  =  0  if  D  =  D. ).  Figure  2  shows  the  wery   small  effect  of  this 
factor  for  the  acceleration  only  case.  Therefore,  if  this  case  is 
a  reasonable  approximation  to  Eq.  (11),  it  is  also  reasonable  to 
write 

K2  x(l  +  H-j)  =  a,  so  that 


t  = 
and  t 


r2¥l  +  Hi>J 

]_2k2(l  +  H^  J 


1/2 


=  1 


(17) 


as  seen  in  Fig.  2. 

b.  Fluid  friction  only:  In  this  case,  Eq.  (11)  becomes 

:2 


and 


t  =  £    TT1 


'k3H3 


1/2 


k3  *  ^H3  +  H4X^  =  a 
\3/2 


so  that 


1/2 


k3H3 


f        H4 


<-k 


-   1 


(Appendix  B) 


3/2 


-  1 


(18) 


Since  H,  is  essentially  equal  to  f  ,  the  friction  coefficient  for 
the  chamber,  it  is  not  likely  to  have  a  value  larger  than  about  0.05. 
I-L  is  of  the  order  of  magnitude  of  1  (in  the  experimental  apparatus 
it  is  about  3.3)  so  that  it  is  unlikely  that  H^/FL  would  be  larger 
than  0.05  in  which  case  tf  (a/k3H3)  ^  =  1.0124.  Thus,  it  is  seen 
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again  that  the  nonlinearity  contributed  by  the  const. (x)  term,  as 

in  the  acceleration  case,  is  relatively  insignificant.  With  this 

approximation,  we  may  write 

2 
k3H3  x  =  a,  so  that 


=  x/4^3 


and  t 


rf 


¥3 


(19) 


This  result  is  all  the  more  useful  in  comparing  the  linearized 
versions  of  the  acceleration  only  and  fluid  friction  only  cases 
Thus, 


1/2 


£  _  ram  time,  acceleration  only 
ram  time,  fluid  friction  only 


linearized 


2k2(l  +  H^ 


k3H3 


(20) 


which  is  a  valuable  ratio  in  analyzing  the  extent  to  which  a  system 
design  is  limited  by  inertia  or  by  fluid  friction  effects. 

c.     Completely  linearized  system:     If  the  fluid  friction  is  com- 
pletely neglected,  along  with  H^,  the  result  is  the  familiar  linear 
damped  spring-mass  expression: 

k2  x   (1  +  H^  +  (k4  +  k?)  x  +  (k5  +  kg)  x  =  a 
with  the  solution 


where 


k5  +  k8 


2       k5  +  k 


1   + 


-CO)     t 

.       n 


l  -  c 


sin     [con/l   -  e2  t  +  a]  J 


(21) 


a  =  tan 


-1 


r 


Vi  - 


-e 


k     +   k 
K4       K7 


wn     =  k2   (1   +  H^  and  2^n  =  k^T^) 
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Since  no  accurate  values  for  kr,  kg,  k,,  and  k-,  were  available  at 
the  time  of  this  writing  (in  the  experimental  apparatus,  these  values 
are  essentially  zero),  this  solution  has  not  been  analyzed  in  any 
great  detail.  In  fact,  our  results  indicate  that  these  damping  and 
spring  terms  do  not  significantly  affect  the  performance  of  an  LPG 
loading  system  when  compared  with  the  similar  effects  present  in  the 
compressible  (gas)  side  of  the  system. 

In  the  following  section  of  this  report  these  closed-form  solu- 
tions are  compared  with  the  experimental  results.  As  has  been  men- 
tioned, the  complete  solution  to  the  governing  differential  equations 
has  been  obtained  by  analog  computer  methods  and  is  discussed  in 
Section  IV  of  this  report. 
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III.   EXPERIMENTAL  PROGRAM 

A.  DESIGN  OF  THE  EXPERIMENTAL  APPARATUS 

After  a  review  of  the  available  background  literature  and  current 
research  efforts,  work  was  begun  to  design  and  fabricate  a  laboratory 
model  capable  of  simulating  the  dynamic  behavior  of  a  liquid  propel  - 
lant  feed  system  for  a  20mm  gun  firing  at  a  rate  of  up  to  1,000 
rounds  per  minute.  The  selection  of  a  20mm  bore  diameter  was  dic- 
tated by  considerations  of  laboratory  capabilities  and  was  not  due 
to  a  desire  to  accurately  simulate  an  operational  system.  The  concept 
criteria  coupled  with  the  unique  requirements  for  an  accurate  labora- 
tory simulation  resulted  in  the  establishment  of  the  following  initial 
design  parameters: 

chamber  volume  (variable)  50-70cc 

maximum  operating  frequency  17hz 

water  flow  rate  (maximum)  75.5gpm 

maximum  chamber  pressure  1 500psi 

desired  maximum  time  to  fill  chamber    15ms 

Due  to  the  time  constraints  involved,  it  was  felt  that  off-the-shelf 
items  should  be  used  whenever  possible.  The  20mm  bore-diameter  chamber 
was  loaded  with  a  brass  slug  weighing  93  grams  which  rode  on  two  graph- 
ite filled  Teflon  sealing  rings.  The  brass  slug,  which  simulated  the 
projectile,  was  cycled  from  the  breech  end  of  the  chamber  to  the  barrel 
end  and  returned  to  the  breech  end,  completing  one  hypothetical  firing 
cycle.  To  facilitate  visual  inspection  of  the  bore  and  slug,  the  test 
chamber  itself  was  fabricated  from  a  three-inch  O.D.  Lucite  cylinder, 
18  inches  long,  bored  to  a  20mm  inside  diameter  and  fitted  with  alumi- 
num end  caps.  The  entire  chamber  assembly  is  held  together  with  four 
tie  rods  and  cap  screws,  as  shown  in  Figs.  3  and  4. 
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Figure  3.    Experimental  Apparatus 


Figure  4.    Gas  Supply  System 
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Because  of  the  desire  to  vary  the  charge  to  mass  ratio,  a  variable 
chamber  volume  was  necessary.  To  accomplish  this  with  one  chamber,  a 
volume  control  retaining  rod  was  designed  into  the  system.  This  brass 
rod,  bored  to  allow  gas  to  pass  its  length,  was  threaded  through  a 
plate  which  was  attached  to  the  barrel  end  cap  holding  the  rod  in  the 
chamber.  The  rod,  which  has  a  Teflon  disc  on  the  end,  not  only  estab- 
lished the  volume  of  the  test  chamber,  but  provided  a  buffer  stop  for 
the  slug  at  the  end  of  its  forward  motion.  Another  Teflon  buffer  was 
affixed  to  the  breech  end  cap  to  cushion  the  slug  in  return  motion. 

The  initial  design  concept  started  a  simulated  firing  cycle  with  a 
slug  at  the  breech  end  of  the  empty  chamber,  as  shown  in  Fig.  5.  This 
is  the  ready-to-ram  position.  The  simulated  propellant  was  then  intro- 
duced, ramming  the  slug  to  the  opposite  (barrel)  end  of  the  chamber  as 
the  chamber  is  filled.  This  was  accomplished  by  applying  gas  pressure 
to  a  power  piston  which  drove  an  injection  piston.  The  injector  piston 
forced  the  propellant  past  a  flow  check  valve  and  into  the  chamber. 
This  placed  the  slug  in  the  ready-to-fire  or  in-battery  position.  In 
an  actual  gun,  the  propellant  would  be  ignited  at  this  time  in  the  cycle, 
Due  to  the  laboratory  constraints,  the  water  expulsion  system  was  used. 

With  the  slug  at  the  in-battery  position,  ram  gas  pressure  on  the 
power  piston  was  vented  and  return  gas  pressure  applied  to  both  the 
power  piston  and  the  muzzle  end  of  the  test  chamber  (Fig.  6).  At  the 
same  time,  gas  pressure  was  applied  to  a  pilot-operated  check  valve, 
the  dump  valve,  to  open  it.  With  the  dump  valve  open  and  reverse  flow 
in  the  feed  line  eliminated  by  the  flow  check  valve,  the  gas  pressure 
on  the  muzzle  end  of  the  chamber  forced  the  slug  to  the  breech  end, 
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expelling  the  propellant  through  the  dump  valve.  At  the  same  time, 
with  the  power  piston  and  injector  piston  moving  to  the  rear,  another 
shot  of  fluid  was  drawn  into  the  injector  from  the  supply  line.  To 
repeat  the  ram  stroke,  the  gas  pressure  was  vented  to  atmosphere. 

Since  the  timing  and  synchronization  of  the  gas  and  liquid  flows 
were  critical  to  the  success  of  the  experiment,  a  single  control  mech- 
anism was  desired.  With  this  in  mind,  an  integrated  gas  and  liquid 
system  was  designed.  It  was  initially  thought  that  balanced  spool 
valves  could  be  used  to  control  one  or  both  ends  of  the  system.  How- 
ever, obtaining  off-the-shelf  stock  for  such  specialized  requirements 
(the  necessary  low  response  times  and  high  flow  rates)  presented  sched- 
uling difficulties  which  were  unacceptable.  It  was  discovered  that 
normally  open  and  normally  closed  solenoid  operated  gas  spool  valves 
were  "readily"  available  for  operating  pressures  up  to  and  including 
150psi  and  could  be  obtained  with  modifications  allowing  operation  up 
to  250psi. 

A  system  was  designed  using  one  valve  to  control  the  ram  gas  pres- 
sure and  one  valve  to  control  the  return  gas  pressure.  By  using  two 
valves,  one  normally  open  and  one  normally  closed,  gas  pressure  could 
be  applied  to  either  the  ram  or  the  return  line  with  the  other  line  be- 
ing vented  by  energizing  or  de-energizing  the  appropriate  solenoids  by 
means  of  a  single  relay.  By  connecting  the  pilot  line  of  the  dump  valve 
to  the  return  gas  line,  the  entire  feed  system  could  be  controlled  with 
one  electrical  switch,  or  each  valve  could  be  energized  and  de-energized 
separately.  This  control  system  would  also  allow  maximum  flexibility 
for  future  efforts  with  little  modification.  Through  an  easily  constructed 
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electrical  timing  circuit,  automatic  ram  and  return  cycling  could  be 
accomplished  for  simulation  of  burst  firing. 

B.  CONDUCT  OF  THE  EXPERIMENT 

The  assembled  experimental  LPG  simulator  is  shown  in  Fig.  7.  A 
linearly  variable  differential  transformer  (LVDT)  was  manufactured  and 
mounted  next  to  the  injector  piston.  This  LVDT  was  attached  to  the 
connecting  rod,  between  the  power  piston  and  the  injector  piston. 
The  volume  of  the  liquid  being  placed  under  pressure  during  each  shot 
was  measured  by  filling  the  system  in  the  ready-to-ram  position  and 
then  draining  it  into  a  graduated  beaker.  By  measuring  the  displace- 
ment of  the  injector  piston  head  with  the  LVDT,  the  volumetric  rate  of 
fluid  injection  during  the  ram  stroke  was  obtained. 

Five  pressure  taps  were  drilled  in  the  test  chamber.  Two  of  these 
taps  were  located  as  close  to  the  breech  end  as  possible,  one  at  20 
degrees  from  top  center  and  the  other  at  20  degrees  from  bottom  center 
in  a  counterclockwise  direction,  as  seen  from  the  breech  end  of  the 
chamber.  The  remaining  three  taps  were  placed  two  inches,  three  and 
one-half  inches,  and  five  inches  from  the  breech  end,  in  line  with  the 
first  bottom  tap. 

It  was  discovered  during  operation  that  when  the  slug  was  returned 
to  the  ready-to-ram  position,  it  would  wipe  some  of  the  liquid  out  of 
these  last  three  taps.   In  the  next  ram  cycle,  bubbles  of  gas  which  had 
replaced  the  liquid  would  distort  the  pressure  response  of  these  taps. 
To  expedite  the  investigation,  these  three  taps  were  plugged  and  only 
the  taps  closest  to  the  breech  were  used. 
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Figure  7.    Chamber  Disassembled 


Figure  8.    Chamber  Assembled 
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A  4-channel  Hewlett-Packard  3960  Magnetic  Tape  Recorder  was  used 
to  FM  record  desired  data  during  system  operation.  Kaman  diaphragm 
type  (lOOOpsi)  pressure  transducers  were  connected  to  the  two  breech 
pressure  taps  on  the  test  chamber.  The  pressure  signals  were  processed 
with  a  Kaman  Digi-Vit  Readout  Unit  which  also  provided  a  visual 
(digital)  display.  Since  the  maximum  output  of  the  Digi-Vit  unit 
was  one  volt,  the  breech  pressure  signals  were  amplified  through 
Sanborn  8875A  Differential  Amplifiers  prior  to  being  recorded.  The 
LVDT  excitation  signal  was  obtained  by  using  a  Hewlett-Packard  Model 
200CD  Oscillator.  The  output  of  the  LVDT  was  passed  through  a  voltage 
divider  circuit  and  then  recorded  on  the  tape. 

A  Brush  Recorder  (Mark  280)  was  used  to  obtain  a  visual  display 
of  the  recorded  data.  By  transcribing  the  desired  signals  on  the 
Magnetic  Tape  Recorder  at  a  tape  speed  of  15  feet  per  minute  and 
playing  them  back  into  the  Brush  Recorder  at  3-3/4  feet  per  minute, 
the  time  scale  of  the  output  was  expanded  by  a  factor  of  four  on  the 
Brush  recordings  (viz.,  from  a  real-time  maximum  of  200mm/sec  to  a 
delayed  time  maximum  of  800mm/ sec. ).  Two  pressures  and  the  LVDT 
reading  were  recorded  for  each  shot. 

Before  each  shot  at  a  new  ram  gas  pressure,  calibration  points  were 
marked  on  the  Brush  chart  to  establish  the  full-scale  deflection  pres- 
sures for  each  transducer  during  the  shot.  The  LVDT  full-scale  deflec- 
tion was  maintained  at  40  divisions  by  adjusting  the  amplitude  and 
frequency  of  the  oscillator.  A  slug  stroke  of  five  inches  was  held 
constant  throughout  the  experiment. 
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To  eliminate  the  possibility  of  particulate  contamination  of  the 
water,  a  25  gallon  tank  of  distilled  water  was  connected  to  the  dis- 
charge of  the  dump  valve.  This  return  line  was  elevated  to  place  a 
hydrostatic  head  of  approximately  one  psi  on  the  discharge  side  of  the 
dump  valve  to  prevent  gravity  draining  and  to  eliminate  air  ingestion 
in  the  chamber  and  connecting  line  when  the  dump  valve  was  open. 

Data  were  obtained  for  nominal  supply  pressures  ranging  from 
50psig  to  the  system  limit  of  220psig  as  indicated  at  the  gas  supply 
system  regulator  (Fig.  8).  Further  detailed  descriptions  of  the  ex- 
perimental apparatus  and  procedures  are  documented  in  a  previous 
report  [22]. 

C.   EXPERIMENTAL  RESULTS 

1.  Chamber  pressure:  Ram  displacement,  chamber  pressure,  and  ram 
gas  pressure  (driving  pressure)  were  recorded  for  driving  pressures 
between  50psig  and  220psig  in  lOpsi  increments.  Typical  results  are 
presented  in  Figs.  9-12.   In  these  figures,  the  recorder  outputs  have 
been  arranged  so  that  a  vertical  line  crosses  each  trace  at  a  corres- 
ponding time  during  the  event.  Event  time  increases  from  left  to  right 
with  one  centimeter  equating  to  12.5  milliseconds.  Ram  time  has  been 
defined  as  the  time  interval  which  starts  when  the  ram  gas  pressure 
trace  deviated  from  the  horizontal  to  the  point  in  time  when  any  dis- 
cernable  volumetric  increase  ceases,  as  shown  in  Fig.  9. 

After  application  of  driving  pressure  to  the  power  piston,  a  time 
lag  was  evident  before  the  piston  was  set  in  motion.  This  time  lag  was 
largely  due  to  static  friction  forces  acting  on  the  pistons  and  the  slug 
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After  an  initial  rapid  increase  of  the  driving  pressure,  these  fric- 
tional  forces  were  overcome,  and  injection  of  the  fluid  and,  hence, 
the  ramming  of  the  slug  proceeded.  The  relatively  gradual  increase 
in  driving  pressure  after  piston  motion  started  has  been  attributed 
to  losses  in  the  upstream  valving  and  connecting  lines  and  the  resultant 
inability  of  the  supply  system  to  provide  an  adequate  flow  of  gas  to  the 
power  piston  chamber. 

During  the  ram  stroke,  the  chamber  pressure  exhibited  some  charac- 
teristics which  were  not  completely  predictable  by  the  analytical  model 
in  its  present  form.  After  an  initial  increase,  the  chamber  pressure 
tended  to  decrease  as  the  slug  preceded  the  fluid  column  down  the 
chamber.  To  facilitate  the  discussion  of  these  pressure  fluctuations, 
one  must  first  observe  the  behavior  of  the  displacement  trace.   In  ad- 
dition to  being  a  measure  of  the  volume  of  fluid  injected,  it  is  an 
indicator  of  fluid  velocity  and  acceleration  behavior.  As  motion  starts, 
the  near-zero  slope  of  the  trace  implies  a  low  velocity  which  increases 
as  expected.  The  rapid  initial  rate  at  which  the  slope  increases  im- 
plies a  rapid  increase  in  the  acceleration.  As  the  ram  stroke  proceeds, 
this  rate  of  increase  of  the  acceleration  diminishes  and  the  curve  be- 
comes more  1 inear. 

As  demonstrated  by  the  analysis,  the  pressure  difference  across  the 
system  is  mainly  influenced  by  inertia  forces  and  viscous  forces.  The 
inertia!  forces  induce  a  pressure  reduction  proportional  to  the  accelera- 
tion. Likewise,  the  viscous  forces  endure  a  pressure  reduction  propor- 
tional to  the  square  of  the  velocity.  Following  the  initiation  of 
motion,  as  the  level  of  acceleration  of  the  fluid  and  the  slug  mass 
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decreases  due  to  friction,  the  difference  in  pressure  between  driving 
pressure  and  chamber  pressure  necessary  to  provide  the  acceleration 
level  decreases  with  an  attendant  increase  in  the  chamber  pressure. 
Initially,  fluid  velocities  are  small  so  that  the  pressure  difference 
is  due  mainly  to  this  large  but  decreasing  acceleration.  As  the  accel- 
eration decreases,  the  chamber  pressure  increases  until  the  velocities 
in  the  system  become  significant  from  a  viscous  friction  point  of  view. 
As  velocity  increases,  the  loss  in  pressure  due  to  these  viscous  effects 
overcomes  the  increasing  pressure  due  to  the  decreasing  level  of  accel- 
eration of  the  moving  solid  and  fluid  masses.  The  net  effect  at  this 
point  is  a  loss  in  chamber  pressure. 

In  several  runs  (Figs.  9  and  10,  for  example),  the  acceleration  of 
the  masses  dies  out  after  a  few  milliseconds.  In  these  cases,  the  ex- 
pected continuous  drop  in  chamber  pressure,  due  to  viscous  effects,  is 
observed  after  the  initial  pressure  peak.  In  several  other  runs  (Fig. 
11,  for  example),  the  acceleration  effects  are  present  throughout  the 
run.  In  these  cases  (typical  of  higher  driving  pressures),  the  velocity 
of  the  moving  masses  increases  through  a  large  part  of  the  event.  Accel- 
eration levels  decrease  initially,  as  before,  with  the  resulting  initial 
rise  in  chamber  pressure.  Viscous  losses  are  then  evident  with  the 
corresponding  decrease  in  chamber  pressure. 

In  a  viscous  fluid,  pressure  decreases  are  proportional  to  the  square 
of  the  velocity  while  the  sliding  friction  losses  of  the  solid  mass  are 
more  accurately  modeled  as  being  proportional  to  lesser  power  of  the 
velocity.  For  a  given  change  in  pressure  drop  across  the  system,  there- 
fore, there  is  a  tendency  for  the  slug  to  accelerate  and  decelerate  more 


41 


rapidly  than  the  fluid.   In  the  present  experiment,  the  slug  apparently 
accelerates  at  decreasing  levels  throughout  the  ramming  process;  no 
deceleration  of  the  slug  was  noted.  Thus,  the  second  and  subsequent 
pressure  decreases  in  the  chamber,  when  they  occur,  are  thought  to  be 
due  to  the  inability  of  the  fluid  column  to  keep  page  with  the  slug. 
This  motion  appears  to  be  oscillatory  in  the  P  curves  of  Figs.  9-12. 

All  of  the  effects  discussed  here  are  influencing  the  motion  of  the 
slug  and  fluid  and  these  remarks  are  not  meant  to  create  an  overly  simple 
interpretation  of  the  complex  interactions  involved.  It  is  also  impor- 
tant to  note  that  the  chamber  pressure  is  measured  at  a  fixed  point  at 
the  breech  end  and  that  the  amount  of  fluid  upstream  of  this  point  de- 
creases with  time.  Thus,  the  length  of  the  fluid  column,  upon  which 
inertia!  and  viscous  forces  are  acting,  decreases  with  respect  to  the 
location  of  pressure  measurement.  This  change  in  length  of  the  column 
is  on  the  order  of  30%  in  these  experiments  and  tends  to  decrease  pres- 
sure drops  due  to  inertial  and  frictional  effects  as  measured  at  the 
fixed  pressure  taps.  The  above  discussion  also  assumes  that  the  driving 
pressure  is  relatively  constant  for  that  portion  of  the  event  in  which 
the  forces  mentioned  are  acting.  The  ram  gas  pressure  did  indeed  behave 
approximately  in  this  manner  at  lower  pressures.  As  ram  pressure  was 
increased,  however,  the  response  deviated  more  and  more  from  the  assumed 
step  input. 

2.  Terminal  behavior  and  ullage  effects:  The  oscillatory  be- 
havior of  the  chamber  pressure  after  the  slug  came  to  rest  was  due  to 
the  expansion  and  contraction  of  the  entrained  gases  in  the  fluid  column 
CFigs.  9-12,  for  example).  At  the  moderate  driving  pressures  initially 
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used,  the  traces  resembled  second  order  damped  responses  and  this 

resemblance  prompted  a  further  investigation  of  this  portion  of  the 

ram  cycle  event. 

Viewing  the  dynamic  portion  of  the  system  as  picture  in  the  sketch 

below,  a  force  balance  on  the  ram  and  injector  pistons  at  the  instant 

the  slug  reached  the  end  of  its  stroke  gives  A  (P  -P„)  =  m  x,-  + 
a  3     pv  r  c '    p  i 

kfp*i  +  kspxi  ^see  Eq'  C1^' 
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From  the  definition  of  the  effective  bulk  modulus,  $  ,  for  a 
mixture  of  liquid  and  ullage  gas,  the  flow  due  to  compliance  in  the 
mixture  may  be  written 

Neglecting  the  mechanical  spring  forces,  k  x  ,  and  combining  these 
expressions  to  eliminate  x., 


P  =  P  - 
c    r 


o-  (m  P  +  k-Pr) 
p  c    fp  c 


V6* 


and,  after  Laplace  transformation,  and  solution  with  Pc(o)  =  0  and 
P  (o)  =  P   ,  the  pressure  variation  is  found  to  be 
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P,(t)  -  Pn  e-«nct 

=  1   +1 sin   /-/■     .     2ti     \  (22) 


where     <j>       =  tan       -  /T~      2 

-  c 
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1/2 
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2A 

P  * 


l/3 


and  ¥-.   =  total  fluid  volume  (liquid  and  ullage) 

Equation  (22)  for  a  damped  linear  system  is  analogous  in  form  to  the 
measured  pressure  variations  following  the  ramming  cycle.  The  natural 
frequency  and  damping  coefficient  can  therefore  be  obtained  from  each 
driving  pressure  curve  by  measuring  the  appropriate  amplitudes  and  the 
period  of  oscillation  (see  Appendix  C).  Using  these  experimental  values 
of  03  and  5  ,  values  for  equivalent  bulk  modulus  and  friction  factor  were 
obtained.  The  variation  in  these  calculated  values  were  within  the  limits 
of  uncertainty  predicted  for  most  cases  [22].  The  percentage  of  entrapped 
(ullage)  gas  in  the  fluid  volume  is  then  obtained  by  using  this  calculated 
value  of  6  .  Trapped  gases  were  found  to  be  on  the  order  of  O.l  to  0.6% 
of  the  liquid  volume  during  these  tests.  As  the  ram  gas  pressure  was  in- 
creased for  successive  runs,  the  oscillation  of  the  chamber  pressure  be- 
came more  erratic  as  can  be  seen  by  comparing  Figs.  9  and  10  with  Figs. 
11  and  12.  The  validity  of  the  results  obtained  for  w    ,  ;,  3  ,  and  kf 
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at  these  higher  pressures  is  therefore  suspect,  even  though  the 
general  form  of  the  curve  repeats  at  all  driving  pressures. 

In  some  of  the  events,  the  chamber  pressure  curve  indicates  a 
pressure  below  atmospheric  immediately  after  the  slug  reaches  the 
end  of  its  travel.  These  initial  low-pressure  oscillations  are 
tentatively  attributed  to  slug/fluid  interactions  as  the  slug  comes 
to  rest  against  the  buffer.  Note  that  following  the  stopping  of  the 
slug,  the  initial  pressure  change  is  always  a  steep  increase,  indi- 
cating that  the  fluid  column  is  initially  compressed  against  the 
stationary  slug. 

It  was  hoped  that  a  value  of  equivalent  bulk  modulus  could  be 
obtained  from  the  displacement  trace  to  correlate  the  value  obtained 
using  Eq.  (22).  However,  electrical  noise  in  the  LVDT  circuitry 
produced  a  signal  with  such  little  definition  that  the  extremely 
small  volume  changes  associated  with  trapped  gas  effects  were  not 
distinguishable. 

3.  Ram  time  measurements:  The  data  displays  of  ram-time  versus 
ram  gas  pressure  (Fig.  13)  portray  some  rather  interesting  charac- 
teristics. The  general  trend  and  validity  of  the  curve  is  confirmed 
by  the  fact  that  the  data  presented  were  taken  on  several  occasions 
and  repeatedly  at  selected  conditions.  It  should  also  be  noted  that 
for  driving  pressures  above  the  moderate  range  (p  «  150psig)  the 
decreases  in  ram-time  are  small  for  increases  in  ram  gas  pressure. 
The  scatter  of  the  data  can  be  attributed  to  the  inherent  uncertainty 
encountered  in  any  experimental  data  taking  system,  coupled  with  pos- 
sible erratic  behavior  of  the  pneumatic  equipment.  Any  uncontrolled 
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variations  in  equivalent  bulk  modulus  for  the  different  tests  did  not 
appear  to  significantly  influence  the  ram-time  results.  An  effort  was 
made  throughout  the  experiment  to  keep  the  amount  of  entrapped  gas  small 
and  relatively  constant.  The  data  points  conforming  to  a  consistent 
trend  would  indicate  that  this  was  accomplished. 

4.  Comparison  with  predictions:  The  experiments  discussed  herein 
have  been  conducted  with  the  following  values  of  the  design  parameters: 


Ac 
c/mt 


Vf*'fc 

k  =k 
Kfp  Kfs 


m 

n 

P 
c 

P 
I 

t 
t. 


chamber  cross-sectional  area,  in 

charge-to-mass  ratio 

injector  diameter,  in.  (=0  >  in  this  case) 

connecting  line  diameter,  in. 

chamber  diameter,  in. 

flow  friction  factors 

mechanical  damping  coefficient,  lbsec/in. 

injector  length,  in. 

connecting  line  length,  in. 

2 
mass  of  injector  piston,  lb  sec  /in. 

2 


0.49) 
0.43) 
1.50) 
0.57) 
0.79) 
0.03) 
.001) 
1.60) 
16.0) 

5.2  x  10"3) 

5.3  x  10"4) 


mass  of  slug,  lb  sec  /in. 

number  of  connecting  lines  (1) 

2 
equivalent  pressure  drop  due  to  drag,  lb/in.  (empirical) 

(controlled) 
(variable) 


ram  gas  pressure,  lb/in. 
time,  sec. 


sm 


ram-time,  sec. 

slug  position,  in. 

maximum  (final)  slug  position,  in. 


(variable) 
(variable) 
(5.00) 
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The  sudden  decrease  in  diameter  from  injector  to  connecting  line  gives 
a  minor  loss  coefficient  (see  Eq.  (4))  of 

.   V*(£)4'" 

The  entrance  to  the  chamber  is  constricted  by  an  orifice-type  reduc- 
tion with  contraction  ratio  (C  )  of  0.92.  The  coefficient  for  this 
is  (with  Cd=0.61) 

k^(iM|)4  =  11-7 

and  T  k  .  =  17.2 
^   mi 
i 

These  design  factors  lead  to  the  following  values  of  the  relevant 

parameters  for  Eqs.  (9)  and  (10): 

k,  =  0 

k2  =  5.41  x  10"3/Pr  sec2 

k..  =  7.36  x  10"3/P  sec2 
3  r 

k4  =  7.85  x  10"4/Pr 

k5  =  k8  =  0 

ky  =  1.02  x  10"2/Pr 

H]  =  3.44 
H2  =  0.39 
H3  =  3.21 
H4  =  0.03 

The  value  of  the  effective  back  pressure,  P,  (see  Eq.  (15))  is  obtained 
empirically  as  is  discussed  in  subsequent  paragraphs. 
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a.  Pressure  and  displacement  variations  with  time. 

Figures  9-12  show  the  results  of  the  analog  computer  solu- 
tions to  the  analytical  model  for  the  conditions  under  consideration. 
The  value  of  P.,  the  effective  back  pressure,  was  set  at  46psi  except 
for  the  P  =  50psi  run  (Fig.  9)  where  it  was  necessary  to  set 
P.  =  41psi.  The  logic  of  these  selections  is  discussed  in  a  subse- 
quent section.  In  obtaining  the  fit  of  the  analytical  model  to  the 
data,  as  was  done  in  Figs.  9-12,  it  was  necessary  to  partition  P. 
into  portions  acting  across  the  injector  piston  and  across  the  slug. 
That  is, 

pd  ■  £ +  r +  pb  ■  pdi+pds 

1    c 

where  Prf.=  ^  and  P   =  /  +  P 
1  c 

In  the  four  cases  illustrated,  it  was  found  that  the  static  drag  was 
almost  entirely  due  to  the  injector  piston.  The  semi-empirical  re- 
sults are  as  follows: 


Figure 

P 
r 

Pd 

Pdi 

Pds 

9 

50 

41 

31 

10 

10 

100 

46 

41 

5 

11 

140 

46 

43 

3 

12 

210 

46 

46 

0 

This  information  is  of  obvious  use  in  pointing  out  the  area  to  seek 
reductions  in  P.. 

The  agreement  of  the  analysis  with  the  experimental  measurements 
shown  in  Figs.  9-12  is  remarkably  good  when  it  is  realized  that  (1) 
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the  step  input  assumed  in  the  model  is  only  crudely  approximated  by 
the  experimental  system  and,  (2),  the  model  does  not  include  compres- 
sibility effects.  The  latter  discrepancy  is  believed  to  be  the  reason 
that  the  analytical  model  does  not  accurately  follow  the  observed 
chamber  pressure  oscillations.  These  oscillations,  caused  by  the 
accordion-like  motion  of  the  liquid  between  the  injector  piston  and 
the  slug,  require  the  accounting  for  compressibility  effects  in  the 
liquid.  A  careful  examination  of  these  results  has  shown  that  the 
deviations  from  the  analytical  model  are  also  present  in  the  dis- 
placement trace,  although  this  is  not  obvious  in  the  figures.  The 
effect  of  the  non-constancy  of  driving  pressure,  P  ,  is  partially 
removed  by  the  semi-empiricism  present  in  the  selection  of  P.,  and 
P  ,,  as  discussed  above.  Finally,  it  should  be  noted  that  the  cases 
presented  in  Figs.  9-12  are  randomly  selected  from  a  large  number  of 
runs.  The  entire  body  of  original  data  is  available  upon  request, 
b.  Ram-time  variations  with  driving  pressure. 

Equation  (14)  gives  a  predicted  ram-time  dependency  based 
upon  a  number  of  simplifying  assumptions.  Repeating  the  expression 
here, 

i  ....-i  ....  M 

A 


where 


tr  =   g-,  cosh   exp 
1  -  fVc  \ 

b'    I  pa  ; 


2m_D_  \  V2         H 


r  c  '       /c_ 

m 


\   r/D  \  2  P.  "I 


i)V2 


and    A  = 


k2  (1  +  H^ 
k3  H3 
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With  the  experimental  conditions  previously  described,  we  have 

1/2    and 


1 

.156 

B' 

<Pr  "  Pd> 

A  =  1.017 

with  the  result  that 

4-      256        ,    v 
tf  =  1/?    (msec) 

(P  -  P A)  ' 

The  agreement  of  this  expression  with  the  experimental  data  is  shown 
in  Fig.  14.  To  obtain  this  agreement,  it  has  been  necessary  to  account 
for  the  approximately  8  millisecond  delay  that  is  observed  in  the  ex- 
periments.  In  addition,  it  is  seen  that  a  good  estimate  of  P.,  for 
the  system  used,  is  46psi.  This  remarkable  agreement  indicates  the 
essential  correctness  of  the  model  when  used  within  the  limits  imposed 
by  the  assumptions.  The  main  uncertainty  in  the  model  is  about  ±10% 
in  the  value  of  H^.  This  uncertainty  is  equivalent  to  about  ±2%  in  tf. 

5.  General  observations.  The  analytical  model  is  seen  to  be  well 
supported  by  the  experimental  data.  It  is  therefore  reasonable  to 
utilize  the  analytical  model  to  extrapolate  the  experimental  results 
to  cases  of  special  interest.  In  particular,  the  analysis  can  be  used 
to  indicate  the  areas  of  redesign  likely  to  be  most  profitable  and  to 
obtain  preliminary  predictions  of  LPG  loading  system  performance  under 
a  variety  of  design  conditions  obtained  through  the  specification  of 
values  of  parameters  different  from  those  analyzed  to  date.  These 
parameters  include  caliber  size,  charge-to-mass  ratio,  and  projectile 
mass  (or  propel lant  density).  Results  of  such  parametric  studies  are 
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reported  in  the  following  section. 

Finally,  it  should  be  pointed  out  that,  as  expected,  the  analytical 
model  only  partially  predicts  the  instantaneous  behavior  of  the  liquid 
and  solid  masses  during  loading.   If  knowledge  of  the  mean  behavior 
with  respect  to  time  is  insufficient  information,  it  is  imperative 
that  the  model  be  expanded  to  take  the  driving  system  (compressed  gas) 
into  account  and  to  allow  for  compressibility  effects  in  the  liquid. 
Some  further  comments  on  this  subject  are  offered  at  the  conclusion 
of  the  report. 
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IV.   DESIGN  CONSIDERATIONS 

In  this  section  the  analytical  model  is  utilized  to  demonstrate 
those  design  decisions  that  are  most  influential  in  determining  the 
performance  of  LPG  propellant  loading  systems.  For  the  purposes  of 
the  following  remarks,  "performance"  is  taken  to  mean  ram-time.  Other 
criteria  (such  as  chamber  pressure  levels  and  system  volume)  are  also 
amenable  to  prediction  by  the  analytical  model  but  ram-time  is  con- 
sidered here  due  to  its  wide  acceptance  as  a  quantity  to  be  minimized 
in  any  LPG  system. 
A.   REDESIGN  OF  THE  LPG  SIMULATOR 

Appendix  D  describes  the  analog  computer  program  that  has  been 
developed  to  solve  the  governing  equations  (9)  and  (10).  This  program 
has  been  used  to  predict  the  behavior  of  the  existing  laboratory  system 
and  to  point  the  way  toward  the  most  profitable  areas  for  redesign  em- 
phasis. 

Figures  15a-e  show  the  computer  outputs  for  slug  displacement, 
velocity,  acceleration,  chamber  pressure,  and  slug  pressure,  respec- 
tively, for  a  range  of  reservoir  pressures  from  100  to  250psig.  The 
values  of  the  parameters  are  those  presented  previously  in  Section 
III.C.4  of  this  report.  The  parameters  k~  and  k3  are  mainly  dependent 
upon  operational  specifications  (slug  mass,  caliber,  charge- to-mass 
ratio,  etc.)  and  are  therefore  not  suitable  for  redesign.  The  other 
k.  parameters  are  negligibly  small  in  the  existing  system.  The 
parameters  H,  and  H3,  representing  inertia  and  viscous  effects, 
respectively,  are  representative  of  two  central  areas  for  redesign. 
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With  H1  =  3.44  and  H3  =  3.21,  we  see  from  Eq.  (20)  that  R  =  1.42  and 
the  system  is  therefore  inertia  dominated  from  a  ram-time  point  of 
view.  The  third  parameter  that  might  be  improved  is  the  static  drag 
on  the  system,  represented  by  the  effective  back  pressure  P.. 

Figures  16a-c  show  the  effects  of  the  parameters  H, ,  H3,  and  P. 
at  a  reservoir  pressure  of  P  =  150psig.  Each  of  these  parameters 
have  been  swept  through  four-fold  reductions  from  the  existing  design, 
From  these  figures  it  is  seen  that  ram-time  decreases  of  3.5,  5.5, 
and  2.5  milliseconds  are  predicted  for  four-fold  decreases  in  P., 
H, ,  and  H3,  respectively.  It  is  apparent,  therefore,  that  redesign 
priorities  should  be  inertia  reduction  (H-,),  static  drag  reduction 
(P.),  and  viscous  loss  reduction  (H3),  in  this  order.  Inspection  of 
the  parameters  H,  and  FL  will  show  that  a  reduction  in  connecting 
line  length  (L«)  and/or  an  increased  connecting  line  diameter  (Dj 
will  decrease  both  H,  and  H3.  A  significant  reduction  in  H,  is 
attainable  through  the  decrease  of  piston  mass  (m  )  or  the  increase 
of  injector  area  (this  term  is  already  rather  small,  however,  rela- 
tive to  the  third  term  in  H-,).  H3  can  be  significantly  reduced 
mainly  through  shortening  of  line  length,  increasing  the  line  diame- 
ter, and,  most  importantly,  reduction  of  minor  losses.  In  the 
current  system,  minor  losses  make  up  about  85%  of  H3  (thereby  point- 
ing up  a  discrepancy  in  terminology).  The  effective  back  pressure 
is  apparently  due  mainly  to  drag  at  the  injector  end  of  the  system. 
This  term  is  also  influenced  (due  to  inadequacies  of  modeling  of  the 
input  pressure)  by  a  non-constant  reservoir  pressure. 
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Figures  17a-e  show  predicted  performance  curves  calculated  by 
assuming  a  two-fold  reduction  in  all  of  the  redesign  parameters. 
That  is,  in  these  calculations,  P.  =  23psig,  H,  =  1.72,  and  FL  =  1.60. 
Referring  to  the  P  =  150psig  curve,  we  see  a  ram-time  reduction  of 
about  7.5  milliseconds.  This  appears  to  be  about  one  millesecond 
more  than  that  calculated  by  simply  adding  the  separate  reductions 
observed  in  Fig.  16  (thus,  the  non-linearity  of  the  problem  is  mani- 
fested). Figure  17  shows  that  if  the  two-fold  parameter  reductions 
are  attainable  through  redesign,  a  ram-time  of  15  milliseconds  should 
be  attainable  with  a  reservoir  pressure  of  about  200psig.  Allotting 
25%  of  the  total  firing  cycle  to  liquid  loading,  this  would  permit  a 
rate  of  fire  1000  rounds  per  minute. 

In  concluding  this  section,  it  should  be  noted  that  all  delays  in 
loading,  other  than  those  due  to  actual  motion  of  the  injector,  liquid, 
and  slug, have  been  omitted  from  this  discussion.  The  8  millisecond 
static  delay  observed  in  our  experiments  would,  of  course,  be  unaccept- 
able in  an  operational  system.  With  this  caveat,  however,  redesign 
recommendations,  in  order  of  importance,  are: 

1.  Reduce  connecting  line  length. 

2.  Increase  connecting  line  diameter. 

3.  Decrease  piston  mass  and  static  drag. 

4.  Seek  out  and  reduce  minor  losses. 

B.   PRELIMINARY  DESIGN  OF  A  5-INCH  SYSTEM. 

In  this  section,  the  analytical  model    is  exercised  in  order  to 
demonstrate  its  utility  in  making  first-cut  design  decisions  applicable 
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'  Figure  17.   Predicted  performance,  LPG  redesign, 
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to  liquid  loading  systems  for  guns  of  different  scale.  Before  pro- 
ceeding, it  is  noted,  once  and  for  all,  that  the  validity  of  these 
predictions  is  vitally  dependent  on  the  extent  to  which  the  basic 
analytical  model  approximates  the  system  under  design.  In  particular, 
any  current  injector  design  will  probably  require  a  more  complex 
mathematical  description  than  that  given  by  Eqs.  (1)  and  (2). 

For  the  purposes  of  this  discussion,  we  shall  rely  mainly  upon 
Eq.  (14)  which  describes  the  semi-linear  case  in  which  H./l-L,  H?/(l+H, ) 
and  (kc+kg)  are  all  negligibly  small.  These  are  desirable  design  con- 
straints and  are  relatively  easily  obtained.  In  any  case,  in  sections 
11. A. 2. a.  and  b.,  the  effects  of  departures  from  these  conditions  have 
been  shown  to  be  small.  Repeating  Eq.  (14)  in  a  slightly  reworked  form 


tf- 


k3H3 


1    V2 


[  Cf )  - 


r    _ 


A  cosh"1   exp  (1/A) 


(23) 


k?(l  +  HJ    9 
where   A  =    -  ,  M — —  =  JT/2 
K3  3 

Note  the  physical  interpretation  of  the  parameter  A:  it  is  a  measure 
of  inertia  effects  relative  to  viscous  effects.  The  function  of  A 
in  Eq.  (23)  is  shown  in  Fig.  18. 

It  is  now  necessary  to  postulate  a  few  design  parameters  descrip- 
tive of  a  5-inch  baseline  LPG  propel lant  loading  system.   It  is  to 
be  emphasized  that  this  baseline  design  need  only  be  approximate 
since  it  is  the  function  of  the  model  to  point  up  profitable  depart- 
ures from  the  baseline.  For  this  example,  the  baseline  system  is 
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described  as  follows: 

Chamber  diameter 
Charge-to-mass  ratio 
Propellant  density 
Slug  mass 
Driving  pressure 


m 

s 
D 
D7  =  °'5> 

sm 


Dc  =  5.0  in. 

c/ms  =0.5 

p  =  1.3  x  10"4lb  sec2/in.4 

ms  =  80  lbm  =  0.207  lb  sec2/in 

Pr  =  250  psig 


p^=0.1, 
r 


tt=  5.0, 


sm 


n  =  2 
D 


D 


=  1.0, 


l 


£ k™-  ■ 20 


The  length  of  slug  travel   can  be  calculated  from 

ms(c/ni  ) 


sm 


pA. 


40.6  in. 
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The  pertinent  parameters  are  found  to  be 

k2  =  1.71  x  10"3  sec2 

k3  =  3.47  x  10"3  sec2 

H]  =  .0625  +  .125  +  6.25  =  6.44 

H3  =  .0009  +  23.44  +  2.46  =  25.90 

Note  the  dominance  of  the  connecting  line  contributions,  6.25  and  23.44, 

to  H,  and  H3,  respectively.  From  these  values,  we  calculate 

A  =  0.141  and  R  =  0.53 
so  that  this  design  is  seen  to  be  dominated  by  viscous  effects.  From 

Fig.  18,  f(A)  =  1.1  and  the  ram  time  is  tf  =  0.35  sec.  With  the  criter- 
ion of  25%  of  the  total  cycle  time  available  for  liquid  loading,  this 
baseline  ram-time  will  allow  a  rate-of-fire  of  about  43  rounds-per-minute. 

There  is  a  limitless  number  of  perturbations  on  the  baseline  design 
that  can  be  performed  in  order  to  examine  their  effect  upon  rate-of-fire. 
Judicious  choice  of  these  perturbations  awaits  the  specification  of  more 
exact  performance  and  design  criteria.  As  examples  of  what  can  be  done, 
however,  it  is  noted  that  this  particular  baseline  design  is  severely 
constrained  by  inertia  and  viscous  losses  in  the  connecting  lines  be- 
tween injector  and  breech.  Thus,  if  the  number  of  1-inch  connecting 
lines  is  doubled  (to  n  =  4),  we  calculate  H,  =  3.31,  and  FL  =  8.32,  so 
that  A  =  0.255  and  R  =  0.714.  The  resulting  maximum  rate-of-fire  is 
about  71  rounds-per-minute,  a  65%  improvement.  As  another  example,  if 
the  injector  diameter  is  constrained  so  that  D  /D.  =  1.0,  the  resulting 
rate-of-fire,  for  the  n  =  4  case,  is  about  55  rounds-per-minute,  a  29% 
decrease. 

As  a  final  demonstration  of  the  analytical  capability  thus  far  dev- 
eloped, the  baseline  5-inch  design  has  been  modeled  on  the  analog  compu- 
ter for  various  driving  pressure  levels.  The  results  are  shown  in  Figs. 
19a-e.  cc 
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Figure  19 
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V.   SUMMARY 

A.  CONCLUSIONS 

An  analytical  model  has  been  developed  to  predict  the  time-dependent 
behavior  of  the  moving  liquid  and  solid  masses  during  an  LPG  loading 
cycle.  This  model  has  been  tested  against  experimental  results  obtained 
with  an  LPG  laboratory  simulator  and  has  been  found  to  be  adequate  for 
the  prediction  of  system  performance.  The  model  is  specifically  con- 
strained in  its  use  by  the  assumption  of  an  incompressible  liquid  and 
a  uniform  step  input  in  driving  pressure. 

An  experimental  LPG  simulator  has  been  designed  and  fabricated  to 
assist  in  verifying  analytical  predictions  and  to  identify  design  and 
performance  phenomena  requiring  particular  attention  in  future  LPG 
developments.  Measurements  made  with  the  LPG  simulator,  together  with 
guidance  provided  by  the  analytical  model,  have  provided  valuable  in- 
sights into  the  behavior  to  be  expected  from  a  given  LPG  feed  system 
design.  Specific  unanswered  questions  relate  to  the  effects  of  large 
trapped  gas  volumes  and  the  existence  and  probable  effects  of  cavita- 
tion in  the  liquid  slug  during  the  loading  cycle. 

The  analytical  model  has  been  shown  to  be  useful  in  the  prediction 
of  the  influence  of  design  specifications  such  as  caliber,  projectile 
mass, charge- to-mass  ratio,  propel  1  ant  properties,  and  system  geometry 
and  scale.   It  is  anticipated  that  the  analysis,  further  substantiated 
by  experimental  observation  and  improved  as  required  for  future  appli- 
cation, will  provide  a  valuable  tool  in  the  design  of  prototype  and 
operational  LPG  systems. 
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B.   RECOMMENDATIONS 

1 .  Experimental 

a.  Incorporate  design  improvements  as  suggested  by  the  analyti 
cal  model.  Test  the  modified  LPG  simulator  to  verify  the 
expected  performance  increases. 

b.  Conduct  tests  to  verify  other  parametric  dependencies  pre- 
dicted by  the  model.  Those  most  easily  evaluated  include 
propel lant  properties  (simulate  a  liquid  propel lant), 
charge-to-mass  ratio,  and  slug  mass. 

c.  Investigate  the  effects  of  large  volumes  (up  to  10%)  of 
ullage  gas. 

d.  Determine  the  existence  of  cavitation  and,  if  it  is  found, 
evaluate  the  effect  of  cavitation  upon  system  performance. 

e.  Develop  the  simulator  to  the  point  that  it  can  be  operated 
in  the  fully  automatic  mode.  Conduct  tests  appropriate 
for  the  evaluation  of  the  system  in  this  mode. 

2.  Analytical 

a.  Continue  to  refine  the  model.  The  difficulty  of  including 
compressibility  effects  may  preclude  this  improvement  un- 
less future  requirements  justify  the  necessary  effort. 

b.  Program  the  governing  equations  for  solution  on  the  digital 
computer.  This  will  be  necessary  if  the  model  is  extended. 

c.  Model  upstream  driving  systems  in  order  to  obtain  accurate 
inputs  to  the  existing  model.  This  will  be  necessary  if 
the  analysis  is  to  provide  performance  predictions  for  the 
fully  automatic  mode  of  operation. 

d.  Monitor  LPG  program  developments  with  a  view  to  subjecting 
emerging  designs  to  evaluation  by  the  model. 
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APPENDIX  A 

NOMENCLATURE 

a  parameter  defined  in  conjunction  with  Eq.  (16). 

A  cross-sectional  area 

C  contraction  coefficient 

C  .  discharge  coefficient 

(c/m  )  charge-to-mass  ratio 

D  diameter 

f  fluid  friction  factor 

F  static  drag  force 

H-...H.  parameters  defined  in  conjunction  with  Eqs.  (9)  and  (10) 

k.p  damping  force  coefficient 

k  spring  force  coefficient 

k  minor  loss  coefficient 

m 

k-. . . kg  parameters  defined  in  conjunction  with  Eqs.  (9)  and  (10) 

L  length 

m  mass 

n  number  of  connecting  lines 

P  pressure 

P.  back  pressure  on  slug 

P.  effective  back  pressure  defined  in  conjunction  with  Eq.  (15) 

R  ratio  defined  in  Eq.  (20) 

t  time 

t-r  final  time,  ram-time 
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x       displacement,  position 

x       final  slug  displacement 
sm  3    r 


effective  bulk  modulus 


e 
l,  damping  factor 


% 


liquid  density 

natural  frequency,  radians/second 


Subscripts 

c       refers  to  chamber 

i       refers  to  low  pressure  (liquid)  side  of  injector;  also 
counting  index 

I  refers  to  connecting  line(s) 

p  refers  to  injector  piston 

r  refers  to  high  pressure  (reservoir)  side  of  injector 

s  refers  to  slug 
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APPENDIX  B 

CLOSED  FORM  SOLUTIONS  TO  THE  GOVERNING 
DIFFERENTIAL  EQUATION 

We  seek  solutions  to  the  equation 

k2x(l  +  H,  +  H2  x)  +  k3x  (H3  +  H4x)  +  (k4+k7)x  +  (x5  +  kg)x  =  a    (Bl) 


x(o)  =  x(o)  =  o 


1.  Semi -linear 


H4      H? 
Let  Tj-  <<  1 ,  i+li  <<  1 ,  k5  +  kg  <<  1 ,  so  that 


k2  (1  +  H-j)  x  +  k3H3x2  +  (k4  +  k7)  x  =  a 
Let  z  =  x,  then  k2(l  +  H, )  z  +  k.,H3z  +  (k4  +  k7)  z  =  a,  and 


dz 


dt 
Ml   +  HJ 


a  -   (k4  +  k7)z  -  k3H3z  "2V'       "1 

and,  after  integration  and  application  of  the  limits 


-2  -1 

tanh    ' 


/T 


-2k^H.z  -   (k.  +  k7h 


T'3 


4  T  7 


/q 


-i  -(k.+k7) 


+  i_tanh-l-^4-7 


fq 


-   k2(l   +  H^ 


where     q  =   (k,  +  k7)     +  4ak3FL 
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Upon  rearrangement  and  solution  for  z  =  x,  we  have 


x  =  AB  tanh  B(t  +  t  )  -  C      (13) 

where  the  constants  are  defined  in  the  main  body  of  the  report. 
Expressions  for  x(t)  and  x(t)  are  obtained  through  integration 
and  differentiation  of  this  expression. 
2.  Inertia  effects  only.  Equation  (Bl)  becomes 

k2  x  (1  +  H]  +  H^x)  =  a 

Multiplying  by  2x 


2x  x  = 


2xa 
k2(l   +  H]   +  H2  x) 


(B2) 


or 


so  that 


mC2\         2a  dx 

d(x   }  "     k^     (1   +  H-,   +  Hgx) 


x  = 

Ho 


2a 


H,x        x-ll/2 


k2H2 


In 


(<<i^) 


1/2 


2(1   +  H,)  2 

Let  u  =     [In   (1   +  1   +  H     x)]  so  that  dx  =  — p u  e     du 


(B3) 


Equation  (B3)  becomes 


2(1  +H,)   2 

-  e   du  = 


H, 


k2H2 


V 


k2H2 


or  i    2a 


dt 


e'  dz  =  t 


and,  since  e 


-u 


H, 


,  upon  rearranging, 


1  + 


1  +  H 


t-(i+£nr  x)  _ 


2k2(l  +  H,)  /I  +  H 


1/2 


D(u) 
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where  D(u) 


2  r[ 

-u   / 

=  e    J 


e  dz  is  Dawson's  integral  [21 J. 


3.  Fluid  friction  only.   In  this  case,  Eq.  (Bl )  becomes 

k3x     (H3  +  H4  x)  =  a 

(H3  +  H4x)1/2  dx  =   (^)1/2  dt 


and,  upon  integration  and  evaluation  of  the  limits, 


Note  that  since  we  have  reduced  the  order  of  the  basic  equation,  the 
initial  condition  x(o)  =  o  cannot  be  satisfied  in  this  formulation. 
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APPENDIX  C 
NATURAL  FREQUENCY  AND  DAMPING 


The  natural  frequency,  u  ,  and  the  damping  coefficient,  r,   , 
can  be  calculated  by  measuring  the  quantities  indicated  on  the  sketch 
below: 


C  = 


V        *  2 

1      ln[a/A]         +  ] 


2tt 


T  V  1    -   C! 


where  u     is  in  rad/sec  and  T  is   in  seconds 


A.     Effective  Bulk  Modulus. 


n 


2        h_lh 

mpVt 


p     t  2 
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00 

n 

natural   frequency,  rad/sec 

(variable) 

m 
P 

.   2 
piston  area,   in 

(1.767) 

2 
1 b-sec 
mass  of  piston,  — r- 

(5.2  x  10"3) 

Vt 

3 
total  volume  of  fluid,   in 

(21.389) 

Be 

equivalent  bulk  modulus,   psi 

(variable) 

The  values  in  parentheses  are  for  this  experimental  setup. 

2 


i  =  0.0354  oj 
e        n 


B.  Friction  Factor,  k~ 


%  '■ 

,kfp 

2AP 

V 

Vt 
3emp 

V 

=  anpl 

°n  ? 

kf  =  0.01035  ojnc 


C.  Trapped  Gas  (Ullage). 


6e  h      6g  Vt 


6  equivalent  bulk  modulus 

B»  liquid  bulk  modulus 

6  gas  bulk  modulus 
9 

Vn  total  volume  of  gas  in  fluid 


For  a  perfect  gas       sQ  »  yP 


where  y  ratio  of  specific  heats 

P    absolute  pressure 


79 


For  air  at  moderate  to  low  pressures,  it  can  be  assumed  that  air 
approximates  a  perfect  gas 

rd1r  =  1.4,   Bg-1.4P 


For  water 

&t   *   3  x  105psi 

Ja  =  i.4P  (3x  1q5  -  %  100% 

Vt         \ 3  x  10DBe   / 

As  previously  reported  in  [22],  the  trapped  gas  volume  thus  calculated 
ranges  from  0.1%  to  0.5%  of  the  total  contained  volume  V.. 
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APPENDIX  D. 

SOLUTION  OF  THE  GENERAL 
GOVERNING  EQUATIONS. 

The  governing  equations,  developed  in  Section  II  of  this  report, 
are  repeated  here: 

1  -  Ps  =  k]  +  k2  x(H1  +  H2x)  +  k3x2(H3  +  H4x)  +  k4  x  +  k5x  +  k6   (D-l) 

Ps  =  k2  x  +  k7  x  +  kg  x  +  kg  (D-2) 

The  solution  of  these  equations,  on  the  analog  computer,  is  developed 
in  this  Appendix. 

The  variables  in  Eqs.  (Dl )  and  (D2)  are  non-dimensional  in  all 

respects  except  for  time.  For  instance,  the  dimension  of  x  as  it 

_p 
appears  in  Eq.  (Dl )  is  sec  .   Non-dimensional ization  with  respect 

to  time  is  now  advantageous  for  scaling  purposes. 

Letting  non-dimensional  time  r  =   gt  and  defining  3  such  that 

32k3  =  1  (D-3) 

we  have,  after  substitution,  elimination  of  P  ,  and  rearrangement 


K2   L  .3 


1  -  (k}   +  k6  +  kg)  -  ^  x'  (H1  +  H2  x)  -  x2(H3  +  H4  x)  - 

(D-4) 

3  <k4 +  h]  \|  -  r9  (k5  +  k8)x 


k 

:2 

Note  that  although  the  notation  has  not  been  changed,  all  variables 
in  Eq.  (D-4)  are  completely  non-dimensional.  The  expression  for 
the  acceleration  has  been  written  in  the  form  of  Eq.  (D-4)  in  order 
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to  illustrate  the  scheme  of  programming  for  the  analog  computer. 
An  auxiliary  relation  for  the  slug  pressure  P  is 

kp 

Ps  =  pif  +  k?ex  +  k8  x  +  kg  (D-5) 

The  chamber  pressure,  P  ,  can  be  developed  by  adding  to  P  the 
pressure  drops  due  to  inertia  and  viscous  effects  in  the  chamber. 
The  non-dimensional  result  is 


Pc  =  Ps  +  Tl  +  T2 


k, 


where  L,  =  fg  x  x2  and  T 2  =  -^      (Mx  x        (D-6) 

Equations  (D-4)  through  (D-6)  are  shown  in  schematic  form  in 
Figs.  D-l(a)  and  D-l(b).  Figure  Dl(c)  illustrates  a  modification 
necessary  when  the  design  is  dominated  by  viscous  effects  (H3>10) 
such  as  in  the  5-inch  baseline  example  of  Section  IV.  It  should 
be  remarked  that  the  logic  of  Fig.  D-l  is  based  upon  unit  program- 
ming (reference  voltage  =  1  unit)  and  that  the  schematic  illustrates 
the  manner  in  which  the  equations  were  programmed  for  solution  on 
two  full -expanded  EAI  TR-20  10-volt  analog  computers.  The  primed 
component  numbers  in  Fig.  D-l(b)  refer  to  the  slaved  computer  and 
numbers  without  primes  refer  to  components  of  the  master  unit.  The 
table  following  Fig.  D-l  lists  the  parameters  determining  the 
potentiometer  settings  in  the  analog  program. 
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Pot  Number 

Parameter 

1 

(k4  +  k7)3 

2 

(k5  +  k8)  Ck3/k2) 

3 

H4/10 

4 

H3/10 

5 

H2/10 

6 

H-,/10 

7 

k2/k3 

8 

k3/iok2 

9 

k7B/10 

10 

k8 

n 

k2/k3 

12 

k9 

13 

l  -  (^  +  k6  +  k9) 

15' 

(k2/k3)(c/ms) 

17' 

f 
s 

Table  D-l 


POT  SETTINGS 
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INITIAL  DISTRIBUTION  LIST 


Naval  Ordnance  Test  Station  10 

Attn:  Code  5032E 
Indian  Head,  MD  20640 

Naval  Sea  Systems  Command  5 

Attn:  NSEA  0331  (2) 

CDR  H.  M.  Effron,  NSEA-65311  (1) 

LCDR  P.  G.  Rainey,  NSEA-032  (1) 

MAJ  C.  L.  Carpenter,  NSEA-03M  (1) 
Department  of  the  Navy 
Washington,  D.  C.  20362 

Naval  Air  Systems  Command  1 

Attn:  Mr.  F.  Marguardt,  NAIR-350B 
Department  of  the  Navy 
Washington,  D.  C.  20362 

Naval  Weapons  Center  3 

Attn:  Mr.  E.  Romero,  4022  (2) 

Mr.  G.  E.  Kovalenko,  3053  (1) 
China  Lake,  CA  93555 

Defense  Advanced  Research  Projects  Agency  1 

Attn:  Mr.  Charles  Lehner  (TT0) 
1400  Wilson  Blvd. 
Rosslyn,  VA  22209 

Ballistics  Research  Laboratories  1 

Aberdeen  Research  and  Development  Center 
Attn:  Mr.  Richard  Comer  (AMXIBL) 
Aberdeen  Proving  Ground,  MD  21005 

Chemical  Propulsion  Information  Agency  1 

Applied  Physics  Laboratory 
Attn:  Thomas  W.  Christian 
8621  Georgia  Avenue 
Silver  Spring,  MD  20910 

University  of  Utah  1 

Attn:  Dr.  W.  A.  Guillory,  Department  of  Chemistry 
Salt  Lake  City,  UT  84112 

Naval  Underwater  Systems  Center  2 

Attn:  Robert  Heffeman  (1) 
G.  E.  Kovalenko  (i ) 
Newport,  RI  02840 
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R&D  Associates  1 

Attn:   Dr.   Ray  Edelman 

P.O.   Box  3580 

Santa  Monica,  CA  90403 

University  of  California  1 

Los  Alamos  Scientific  Laboratory 

Attn:  Dr.  Dan  Butler  (Code  T-3) 

P.O.  Box  1663 

Los  Alamos,  NM  87554 

Pulsepower  Systems,  Inc.  1 

Attn:  Mr.  L.  Elmore 
815  American  St. 
San  Carlos,  CA  93555 

General  Electric  Co.  1 

Armament  Systems  Dept. 
Attn:  E.  Ashley 
Burlington,  VT  05401 

General   Electric  Co.  1 

Attn:  Mr.  J.   Has  kins   (Rm.   1046) 
100  Plastics  Ave. 
Pittsfield,  MA  01201 

Stanford  Research  Institute  1 

Attn:  Mr.  G.  A.  Branch 
333  Ravenswood  Ave. 
Menlo  Park,  CA  94025 

Library  2 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Dean  of  Research  Administration  2 

Naval  Postgraduate  School 
Monterey,  California  93940 

Director  *2 

Defense  Documentation  Center 
5010  Duke  Street 
Alexandria,  Virginia  22314 

Chairman,  Department  of  Mechanical  Engineering  1 

Naval  Postgraduate  School 
Monterey,  California  93940 

R.  H.  Nunn,  Associate  Professor  of  Mechanical  Engineering       10 
Naval  Postgraduate  School 
Monterey,  California  93940 
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